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@ Tracking system to fdllow the position and orientation of a device with radioftaquency fields. 



(g) A tracking system in which radiofrequency 
signals emitted by an invasive device (120) such 
as a catheter, are detected and used to measure 
the position and orientation of the invasive 
device (120). The invasive device (120) has a 
transmit coQ attached near its end and is driven 
by a low power RF source to produce a dipde 
electromagnetic field that can be detected by an 
array of receive coils (160) distributed around a 
region of interest The position and orientation 
of the device as detemnined by the tracking 
system are superimposed upon independently 
acquired Medical Diagnostic images, thereby 
minimizing the radiographic exposure times. 
One or more invasive devices (120) can be 
simuKaneousty tracked. 
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BACKGROUND OF THE INVENTION 

The present invention relates to medical procedures in which an invasive device such as a catheter, guide 
wire, biopsy needle, endoscope, laparoscope or the like is inserted into a body, arid more particularly concerns 
5 the tracking of such a device without the use of X-rays. 

X-ray fluoroscopes are used routinely to monitor the placement of Invasive devices during diagnostic and 
therapeutic medical procedures. Conventional X-ray fluoroscopes are designed to minimize X-ray dosage. 
Nevertheless, some procedures can be very long and the accumulated X-ray dose to the patient can become 
significant. The long term exposure of the attending medk:al staff Is of even greater concern since they par- 
10 ticlpate in these procedures regularly. Consequently, It is desirable to reduce the X-ray dose during these pro- 
cedures. 

Another limitation on the use of X-ray fluoroscopes is that the technk^ue Is projective In nature and pro- 
duces a single two-dimensional Image, information concerning the depth of an object within the fleld-of-vlew 
is not available to the operator, it is often desirable to obtain this Information during Invasive procedures. 

16 

SUMMARY OF THE INVENTION 

The present invention provides a tracking system for use in an imaging system for following the iocatton 
of at least one invasive device within a subject, comprising: 
20 (a) means for creating an electromagnetic field of known geometry within the subject; 

(b) an RF receiver for detecting the electromagnetic field at a plurality of M selected locattons; and 

(c) a tracking means for computing position and orientation of the electromagnetic field creation means 
at the M selected locations, responsive to the detected electromagnetic field. 

Tracking of catheters and other invasive devices without X-rays is accomplished using RF transmitters 
25 and receivers. An invasive device such as a guide wire, catheter, endoscope, laparoscope or btopsy needle is 
modified by attaching a small transmit coil near Its end. This transmit coH is driven by a low power RF source 
and creates a dipole electromagnetic field. This dipoie field induces currents and voltages In an array of receh/e 
coils distributed around a region of interest These voltage signals from the receive coils are digitized and sent 
to a tracking computer for analysis. The tracking computer utiib^es non-linear iterative methods to solve for the 
30 posltton and orientation of the transmit coil. This positional information Is then superimposed on a video image 
of the region of interest. Simultaneous detection of multiple transmit coils to better characterize the position 
and orientation of the Invasive device can be accomplished using both time and frequency multiplexing 
schemes. 

Features of the present Invention are to provide a system for tracking an Invasive device without using X- 
36 rays; to provide a method of tracking an Invasive device which minimizes any X-ray does to the patient and 
medical staff; and to provide an interactive Image of an Invasive device superimposed upon another medical 
diagnostic Image. 

BRIEF DESCRIPTION OF THE DRAWINGS 

40 

The features of the invention believed to be novel are set forth with particularity in the appended claims. 
The invention itself, however, both as to organization and method of operation, together with further objects 
and advantages t hereof, may best be understood by reference to the following description taken in conjunctksn 
with the accompanying drawing in which: 
45 Figure 1 isaperspectiveviewofoneembodimentofthepresentlnventioninoperationtrackingthelocatton 
and orientation of an invasive device in a patient 

Figures 2A and 2B together are a schematic block diagram of a radlofrequency tracking system according 
to the present Invention. 

Figure 2C is a diagram Ulusb^ting how Figures 2A and 2B are assembled. 
so Figure 3 is a vector representation of an electromagnetic dipole located at the origin. 

Figure 4 Is a vector representation of an electromagnetic dipole located at a position In space other than 
the origin. 

Figure 5 Is a flowchart of the method used to compute the position and orientation of an invasive device. 

65 DETAILED DESCRIPTION OF THE INVENTION 

In Figure 1, a support arm 101 capable of being rotating about at least one axis 102 and translated by 
gantry control nrieans 70 is shown. Support arm 101 holds an X-ray source 103 that emits a substantially ool- 
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limated beam of X-rays 104 suitable for X-ray imaging and X-ray fluoroscopy. Support arm 101 also holds an 
X-ray detection means 1 05 aligned with the propagation direction of X-rays 104 emitted by X-ray source 103. 
X-rays 104 penetrate a subject support table 110 and a subject 112. An invasive device 120 is placed inside 
the subject by an operator 140. The location of the invasive device 120 is visible on the display of an X-ray 

5 image on a display monitor 1 51 of display means 150 driven by a tracking/display unit 108. In fluoroscopic us- 
age, this image is acquired and displayed several (12 to 60) times a second. 

According to the invention, a plurality of M RF receive colls 160 are placed about the subject, In the pre- 
ferred embodiment, receive coils 160 are attached to X-ray detection means 105. Invasive device 120 Is modi- 
fied to incorporate a small RF transmit coil (not shown in Figure 1 ). The transmit coils may be attached to several 

10 Invasive devices, with at least one coli per invasive device to determine its position and at least two coils per 
Invasive device to determine its orientation. Tracking/display unit 108 provides power to the RF transmit coll 
to create a dipole electromagnetic field which is detected by RF receive coils 160. The signals detected by re- 
ceive coils 160 are used by tracking/display unit 108 to calculate the position and orientation of the b-ansmit 
coil (and therefore invasive device 120). The calculated positbn of Invasive device 120 is displayed by super- 

16 position of a symbol 1 52 on an X-ray image appearing on video monitor 1 51 . 

Following the preferred procedure, operator 1 40 initiates acquisition of X-ray image only when it is deemed 
necessary, in order to minimize X-ray dose to subject 112 and operator 140. The instantaneous tocatton of in- 
vasive device 120 is updated several times per second (ideally 12 to 60 times per second) and provides an 
approximation of the fluoroscopic image of invasive device 120 that operator 140 would expect to see with a 

20 conventional X-ray fluoroscopic system. 

Tracking/display unit 108 is comprised of an RF transmitter 5 and an RF receiver 7 as shown in Fig. 2A, 
and analog-to-digital (A/0) converters 48a, 48b, 48m. a tracking computer 50, a frame grabber 54, and a su- 
perposition means 56 as shown in Fig. 2B. RF transmitter 5 employs a master oscillator 10 that generates a 
signal at a selected frequency. This signal is propagated to a plurality of N frequency offset means 20a, 20b, 

25 20n which generate a plurality N signals of selected different firequencies. Each frequency offset means 20a, 
20b, 20n propagates its signal to a gating means 21a, 21b, 21n, respectively, which either passes the signal 
to an amplifier means 23a, 23b, 23n. respectively, or blocks the signal. Amplifier means 23a, 23b, 23n boost 
the signals by a selected gain G1 and drive transmit coils 30a, 30b, 30n, respectively. In the preferred em- 
bodiment, a number N of these transmit coils 30a, 30b, 30n are situated on invasive device 120. 

30 The signals generated by the transmit coils are detected by a plurality of M receive coils 40a, 40b, 40m 
placed at known locations and with known orientation about the subject Each receive coil 40a. 40b, 40m de- 
tects signals emitted by all transmit coils. The amplitudes and phases of these detected signals are a functton 
of the relative placement and orientattons of the transmit and receive colls. The signals detected by each re- 
ceive coil 40a, 40b, 40m are propagated to low-noise amplifiers 42a, 42b, 42m, respectively, in RF receiver 7, 

35 where they are amplified by a selected gain factor G2. The amplified signals are passed from low-noise am- 
pllf iera 42a, 42b, 42m to quadrature phase detectors 44a, 44b, 44m, respectively, where they are mbced with 
a reference signal from a receive frequency offset means 43 that is driven by master oscillator 10. Mbdng two 
signals In each quadrature phase detector results In a signal having a component at a frequency equal to the 
sum of the Input frequencies, and a component at a frequency equal to the difference of the input frequencies. 

40 The component of interest in the preferred embodiment of this invention is the component equal to the differ- 
ence of the input frequencies. The signals are propagated in quadrature fashion (l.e. as pairs of signals having 
a 90 degree phase difference) to filters 46a, 46b, 46m, respectively, where the low fiBquency component Is 
selected and propagated to A/D converters 48a, 48b, 48m, respectively. The A/D converters convert the low 
frequency signals in each quadrature pair to digital form. This digital informatton is sent to tracking computer 

45 50 through a data bus 51. The tracking computer calculates the positions and orientations of the N b^nsmit 
coils using the digitized signals derived from the M receive colls. The calculated positions and orientattons of 
the N transmit coHs are transformed to display coordinates. 

X-ray imaging and fluoroscopy system 52 is used to generate a video signal that Is propagated to frame 
grabber means 54 which captures a single image from X-ray system 52. Frame grabber means 54 propagates 

50 the single X-ray image in video form to superposition means 56 which overlays a symbol 152 on the video 
signal supplied by frame grabber means 54. The composite video signal is propagated to a suitable display 
means 1 50 such as video monitor 1 51 shown in Figure 1 . The X-ray image will detect the invasive device. The 
tracking computer 50 is initialized by placing the invasive device at an origin marked on the table 110, and 
setting the position to zero. The X-ray system, likewise is adjusted to coincide with the origin marked on the 

55 table. 

Tracking computer 50 communicates with a control computer 60 (Fig. 2A) through an interface connection 
59. Control computer 60 is also Interfaced through a control bus 62 to transmit firequency offset means 20a. 
20b, 20n. gating means 21a, 21b, 21n, transmitter amplifier means 23a, 23b. 23n, receive fi«quency offset 
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means 43, and filters 46a. 46b, 46m. Furthermore, tracking computer 50 is interfaced through an interface 
connection 75 to a gantry control means 70 which is capable of changing the relative position and orientation 
of the subject and the x-ray detection means 105 (Fig. 1). Control computer 60 is responsive to a timing signal 
firom master oscillator 10. 

5 In the preferred embodiment of this invention transmit coils 30a, 30b, 30n are placed on invasive device 

120 illustrated in Figure 1. Reciprocity between receive and transmit coils exists, and placement of receive coils 
40a, 40b, 40m on Invasive device 120 and placement of transmit coils 30a, 30b, 30n outside the subject Is pos- 
sible. 

In the illustrated embodiment of the Invention, a minimum of N^i transmit colls and M^S receive coils is 

10 required to unambiguously determine the location and orientation of invasive device 120. It can be advanta- 
geous, however, to have N>1 transmit coils to provide location and orientation for multiple points on the invasive 
device and/or multiple invasive devices. 

Several methods for detecting signals from multiple transmit coils are possible. One method requires that 
only one of the N gating means be allowed to propagate signal at any Instant. Selection of the propagating 

IS gating means Is made by control computer 60. Control computer 60 notifies tracking computer 50 of a gating 
means selected. The process is repeated for each of the N coils. Thus, tracking computer 60 is able to calculate 
the position of N transmit coils. 

An alternative embodiment requires that all N transmit coils be active at once, each transmitting at a dif- 
ferent frequency. If all N transmit firequencies are within the selected bandwidth of each filter means, then a 

20 collection of L data points can be acquired from each AID converter means. The data points are denmjltiplexed 
by Fourier or Hadamard b^nsformation to separate the individual frequency components arising from each 
transmit coil. Alternatively, M receivers can be constructed for the N transmitters if each transmit frequency is 
within a bandwidth of the filter means of the M receivers. 

Turning now to Figure 3 a vector representation of an electromagnetic dipole 200 located in a three- 

2S dimensional coordinate system comprising an X axis, Y axis and Z axis, and having an origin 201. is shown. 
The strength of the electromagnetic field generated by the dipole at a given positton 205 (given as x, y and z) 
in three-dimensional space is a function of position 205, the orientation of the dipole, here defined by rotation 
angles 9 and ^ and the physical constant Mo known as the permeability of free space, and can be expressed 
as: 

30 

B(w)«^Hi^(M,(2xl.z2-y2) + 3Myxy +3M,xzyi 

^ +(My(2yi-x2.i2) + 3M»xy +3M,yz)j 

+ (Mx(2z2-x2.y2) + 3M»xz +3Myzy)k] 
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In this equation the electromagnetic field at a selected positton 205 in space is divkied into three ortho- 
gonal components defined by the unit vector quantities i. j and k. R represents the distance between the lo- 
catton of the dipole and the selected position, an d is defined as: 

R = Vx2 + y2 + z2 [2] 

The quantities M„ My and Mz represent the vector components of the unit dipole along the x, y and z axes 
and can be expressed as: 

Mz = cos (9) [3] 
My » sin (<|^) sin (9) [4] 
M, " cos (<|)) sin (9) [5] 
where 9 and ^ are the angles shown in Fig. 3. 

In the present lnventk)n it is convenient to translate the location of the dipole to a position other than the 
origin, as illustrated in Figure 4. Since the coordinate system is simply translated and not rotated, the rotatksn 
angles in the new coordinate system 9' and ^' are identical to the original rotation angles 9 and ^. The translated 
origin is 201' (xq, yo. zg). The translated dipole 200' creates an electromagnetic field at a selected receive coil, 
i, In space at position 205 (X|, yi, Z|) which can be calculated from equatton 1 using the following substitutions 
forx, yandz: 

X = X, - xo [61 



4 



EP 0 531 081 A1 



y a yi - yo [7] 
z .= z, - zo [8] 

Each receive coil is positioned at a predetermined location, with receive coi 1 being located at (xi, yi, Zi), 
receive coil 2 at (X2, y2» Z2). etc. Receive coil 1 experiences an electromagnetic field of flux density Bi at location 
5 (Xi, yi, Zt) from the transmit coil transmitting from point (Xq, yo. Zo). translated from the origin. 

The same transmit coi at the same point (xo. yo, Zo). causes receive coil 2 to experience a magnetic field 
of flux density B2 at a location (X2, y2, Z2). This is true for all receive coUs at a given Instant. 

Since x^XpXo . y^rVo * and z^ztZo , x. y and z in equation 1 can be replaced by known quantities resulting 

in: 

10 

Bi(xi,y i.Zi) » ^[((M. {2(xi-xoy^ - (zj-zop - (yi-yo)?| + 3 My(xi.xoKyi-yo) + 3 Mi Xi-xoXzi-zc))! 

16 + (My {2(yi-yo)2 -(xj-xoF -( W) + 3 M. (xi-xo)(yi-yo) + 3 M^iyj-yoM) j 

+ {(Mz{2(zi.zoJ^-(xi.xo)^^yi.yoy^} + 3M«(xi.xo)lz^^^^ ^-3 My(zi.zo)(yi-yo]> k] (^j 

20 

for each of i=1 to M receive coils. These equations can be generalized as: 

B, = f(e. Xo. yo, Zo) 
82 = f(e, (>. Xo. yo. Zo) 
Bs = f(e. Xo, yo. Zo) 

25 B4 » f(e.*.Xo.yo.Zo) 



B„ = f(0. 4>. Xo, yo. Zo) 

30 which is a set of M equations with 5 unknowns. These equations are solvable provWed M s 5, meaning that 
there must be at least 5 receh/e coils. 

Tracking computer 50 (Fig. 2B) Is used to solve the (x. y, z) posltton data for the transmit cols acting as 
a dipole. A flow diagram of a method of solving for the position data is shown in Figure 5. Step 301 serves as 
the entry Into the method. Step 303 Initializes variables to appropriate values. 

35 Step 304 provides an Initial guess for the location and orientation (0, x©. yo, Zo) of the N transmit coils 
being tracked. This initial guess Is obtained by placing the invasive device at a predetermined position which 
is marked on table 110 (Fig. 1) as being (xo.yo.Zo) location, and aligning the invasive device at a predetermined 
orientation (9, ^) also marked on the table at the beginning of the tracking process. 

At step 305. data are acquired from the M receive coils. At step 307, the N components contained in each 

40 of the M sets of data acqured by the receive coiis are separated. This separation can be accomplished by 
Fourier or Hadamard transfbrmatton or any other appropriate methods. At step 309. the position of each of 
the N transmit ooHs Is computed by using equatton [9]. Step 309 can be performed by any suitable mathemat- 
ical algorithm, although our preferred method is an Iterative non-linear optimization such as a multi-dimen- 
sional Newton-Raphson procedure as explained \nNumerk:al Recipes, The Artof Scientific CompuUng,bymu 

45 nam H. Press. Brian P Flannery. Saul A. Teukolsky and William T. Vetterling. Cambridge University Press. 1986. 
chapter 9.4 -Newton-Raphson Method Using Derivatives" hereby incorporated by reference. The Newton- 
Raphson procedure is applied to minimize the difference between the data acquired from the M receive coils 
and the predicted data calculated by using equattons [1-9] for the best guess of the N transmit coi positions 
(xo,yo.zo) and orientations (9.*). At step 311 of the algorithm, the calculated position (xo.yo.Zo) of the N transmit 

50 coils 30a, 30b, 30n is displayed by the superpositton means 56 and display means 150 (Rg. 2B). Step 313 of 
the algorithm determines whether the tracking process is complete. If the tracking process is complete, step 
315 of the algorithm Is taken and the process stops; otherwise, a new guess of the position (xo,yo.Zo) and ori- 
entation (0.<^) for each of the N transmit coils is made at step 317 of the algorithm. The presently preferred 
method to make the guess at step 317 is a linear extrapolatton of the positton and orientation based on the 

55 two Immediately prior positions and orientations calculated for each of the N colls. After step 317 of the algo- 
rithm Is complete, the process continues with acquisitton of new data at step 305. 

Referring now to Rg. 1. the present Invention Is also concerned with the automatic placement and align- 
ment of the subject, by use of a support arm 1 01. within a desired regton around Invasive device 120. This is 
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accomplished in the present instance by transferring the calculated position of the invasive device from track- 
ing computer 50 (Fig. 2B) to a means 70 for controlling the position and orientation of support arm 101 in re- 
lation to support table 110. An X-ray image can be also initiated whenever invasive device 120 enters a region 
of the subject for which an additional X-ray image is required. This embodiment frees the operator from the 

5 task of Iceeping the invasive device 1 20 within the f ield-of-view of the X-ray system 52 and potentially reduces 
the number of assistants that the operator requires. 

Radiologic images upon which the symbol representing the Invasive device is superimposed may be ob- 
tained by means other than X-rays. Data obtained with Magnetic Resonance scanners, Ultrasound scanners, 
Positron Emission Tomography scanners and the like can be used in place of X-ray system 52. 

10 While several presently preferred embodiments of the novel radlofrequency tracking system have been 
described in detaii herein, many modifications and variations will now become apparent to those skilled in the 
art It Is, therefbre, to be understood that the appended claims are intended to cover aO such modtflcatlons 
and variations. 



Claims 

1. A tracking system for use in an imaging system for following the location of at least one invasive device 
within a subject, comprising; 

20 (a) means for creating an electromagnetic field of known geometry within the subject; 

(b) an RF receiver for detecting the electromagnetic field at a plurality of M selected locations; and 

(c) a tracking means for computing positton and orientation of the electromagnetic field aeatton means 
at the M selected locations, responsive to the detected electromagnetic field. 

2. The system of daim 1 further comprising; 

(a) means for acquiring a medical diagnostic image; 

(b) superposition means for superimposing a symbol at the computed position representing the elec- 
tromagnetic field creation means on the medical diagnostic image; and 

(c) display means for displaying the superimposed image on the medical diagnostic image. 

3. The system of claim 1 , wherein the means for creating the electromagnetic field comprises: 

a) a master oscillator for setting a standard timing signal; 

b) a control computer coupled to the master oscillator and responsive to the standard timing signal; 
and 

c) a plurality of N transmit branches, wherein N £ 1, each transmit branch comprises: 

i. transmit frequency offset means for aeating an RF transmit signal of a frequency determined by 
the control computer based upon the standard timing signal, 

ii. a gating means coupled to the transmit frequency offset means and responsive to the control com- 
puter for passing or not passing the RF transmit signal, 

111. an amplifier responsive to the control computer for receiving the RF transmit signal firom the gat- 
ing means and amplifying the signal with a gain determined by the control computer, and 
Iv. a transmit coll coupled to the amplif ierfor creating an electromagnetic field when the RF transmit 
signal Is received. 

4. The system of claim 1 , wherein said means for creating the electromagnetic field of known geometry with- 
in the subject comprises means for generating an electromagnetic dipole. 

5. The system of dalm 3, wherein at least one transmit coll Is attached to one of the group consisting of a 
guide wire, a catheter, an endoscope, a laparoscope, and a biopsy needle: 

^ 6. The system of dalm 3, wherein at least one transmit ooil is incorporated into a surgical devtoe. 

7. The system of dalm 3, wherein at least one transnrUt ooil is attached to a therapeutic device. 

8. The system of dalm 1, wherein the RF receiver comprises: 

55 a] receive frequency offset means for generating a reference frequency signal; and , 

b) a plurality of receive branches, each branch comprising: 
i. a receive coil for sensing the electromagnetic field, 

11. an amplifier coupled to the receive coil for produdng an amplified signal having high and low fre- 
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quency components in response to said receive coil sensing said electromagnetic field, 

iii. a quadrature phase detector coupled to the amplifier and the receive frequency offiset means 
for sampling the amplified signal and comparing it to the reference frequency signal, 

iv. a filter for filtering out the high frequency components and passing the low frequency component 
of the sampled signal, and 

V. an A/D converter coupled to the filter for creating a digital signal representing the electromagnetic 
field within the subject as detected by the receive coll. 



9. The system of daim 8, wherein the receive coO is attached to an invasive device comprising one of the 
10 group consisting of a guide wire, a catheter, an endoscope, a laparoscope and a surgical device, and said 

means for aeating the electromagnetic field of known geometry Is located in predetermined positions 
within the subject. 

10. The system of daim 2, wherein the means for acquiring a medical diagnostic image comprises an X-ray 

system. 

11. The system of daim 1 0, wherein said X-ray system indudes means for providing a selected X-ray image 
from a fluoroscopic sequence. 

12. The system of daim 2. wherein the means for acquiring a medical diagnostic Image comprises a Magnetic 
^ Resonance system. 

13. The system of daim 2, wherein the means for acquiring a medical diagnostic image comprises a Positron 
Emission Tomography system. 

25 14, The system of daim 3. induding means for permitting only a selected one of the transmit branches to cre- 
ate an electromagnetic field at any given instant in time. 

15. The system of daim 3. induding means for permitting the N transmit branches to simultaneously create 
N electromagnetic fields for detection by the M receive branches, said tracking means being responsive 
30 to said N elecferonrmgnetic fields fbr computing N positions and orientations. 

18. The system of daim 8. induding means for permitting the M receive branches to simultaneously detect 
N electromagnetic fields from the N transmit branches, said tracking means being responsive to said N 
electromagnetic fields for computing N positions and orientations. 

^ 17. The system of daim 16. wherein the tracking means incorporates a demultiplexing algorithm comprising 
one of the group consisting of a Fourier transform and a Hadamard transforn). 

18. The system of daim 3, including means for permitting a plurality of the transmit branches to create elec- 
tromagnetic fields of different RF frequendes at a given instant, and a plurality of detection means for 

^ detecting said electromagnetic fields of different RF frequendes, each of said detection means being re- 
sponsive to different selected radtofrequencles, respectively. 

19. The system of daim 3, wherein a plurality of invasive devices are tracked with each invaswe device having 
at least one transmit coa. 
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